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ABSTRACT 

Aims. To determine the effect of the Hall term in the generalised Ohm's law on the damping and phase mixing of Alfven waves in the 

ion cyclotron range of frequencies in uniform and nonuniform equilibrium plasmas. 

Methods. Wave damping in a uniform plasma is treated analytically, whilst a Lagrangian remap code (Lare2d) is used to study Hall 

effects on damping and phase mixing in the presence of an equilibrium density gradient. 

Results. The magnetic energy associated with an initially Gaussian field perturbation in a uniform resistive plasma is shown to decay 

algebraically at a rate that is unaffected by the Hall term to leading order in k^Sj where k is wavenumber and i5, is ion skin depth. 

A similar algebraic decay law applies to whistler perturbations in the limit k^Sj » 1. In a nonuniform plasma it is found that the 

spatially-integrated damping rate due to phase mixing is lower in Hall MHD than it is in MHD, but the reduction in the damping rate, 

which can be attributed to the effects of wave dispersion, tends to zero in both the weak and strong phase mixing limits. 
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1. Introduction 

The interaction between Alfven waves and plasma inhomo- 
geneities forms a well-studied and important area of research, for 
both laboratory and astrophysical plasmas. One process which 
arises as a result of this interaction is Alfven wave phase-mixing. 
Early studies of Alfven wave phase-mixing demonstrated a po- 
tential for significantly enhanced plasma heating. In particular, 
|Heyvaerts & Priest ( 1983| l proposed the phase-mixing of Alfven 
waves as a potential solar coronal heating mechanism through 
enhanced wave-dissipation. They outlined that for a magneto- 
hydrodynamic (MHD) treatment of initially planar shear- Alfven 
waves, propagating independently on individual magnetic sur- 
faces, large differences in phase are quickly generated between 
waves on neighbouring field fines, as a result of variation in 
Alfven speed across tfie field. These phase-differences generate 
progressively smaller scales, and dramatically enhance the ef- 
fects of viscous and Ohmic dissipation, in the locations where 
the Alfven speed gradient is steepest. 

This initial concept has been subsequently adapted for a va- 
riety of problems using MHD, based on the premise that the 
waves in question aiise as a result of an infinite series of bound- 
ary motions at the photosphere. For example, this treatment has 
been used to investigate the heating of open magnetic field lines 



been studied as a source of non-linear coupling to other wave 
modes (see, e.g. |Nakariakov et al.|1997| [Botha et al.|2000| l. 



Of particular interest for this paper is the work of Hood et al. 
( 2002[ l. They note that an infinite series of boundary motions 
is unrealistic and instead investigate the effect on the amplitude 
damping rate due to phase-mixing when the waves are driven by 
only one or two initial impulsive motions at the boundary. 

Recent studies have also begun to move away from the orig- 
inal MHD treatment, instead focussing on full kinetic descrip- 
tions of a plasma undergoing phase-mixing in a colUsionless 
regime, as a potential mechanism for electron acceleration (see. 



e.g. 'Genot et al.'2004', 'Tsikl auri et al.|2005t [Tsiklauri & Haruki 
2008; Bian & Kontar 2010) . On the assumption that the wave- 
lengths of interest (A) are small compared to the particle mean 
free path A^/p, Tsiklauri et al. ( 2005 | l and Bian & Kontar (2010 1 
model the corona as a coUisionless plasma and cite Landau 
damping as their primary wave dissipation mechanism. Such 
damping is strongly suppressed if /I » A^/p ( Ono & Kulsrud 



1975), as in the case of propagating EUV disturbances with pe- 
riods of tens or hundreds of seconds (see e.g.|De Moortel|2009 



under various conditions, e.g. by Parker ( 1991 1; Hood et al. 
( [TWTI i; [DeMoortel et al.| ( [T999l [2000) 1. Phase-mixing has also 



and references therein), and a fluid model is then appropriate for 
the coronal plasma. On the other hand for waves with frequen- 
cies approaching the ion cyclotron frequency, typical coronal pa- 
rameters correspond to classical (Spitzer) collisional mean free 
paths that exceed A, suggesting instead the vaUdity of the col- 
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lisionless approach. However, |Craig & Litvinenko] ( |2002[ ) have 
shown that it is not appropriate to use classical resistivity under 
flaring conditions because it implies current scale lengths that 
are several orders of magnitude shorter than A^fp, and moreover 
is inconsistent with the electric fields required to account for the 
observed acceleration of protons to tens of MeV on timescales 
of the order of one second (Hamilton et al. 2003[l. Craig & 



Litvinenko proposed that the effective resistivity under flaring 
conditions (specifically in a reconnecting current sheet) is de- 
termined by turbulence arising from electron-ion drift (e.g. ion 
acoustic) instabilities, and deduced that this effective resistivity 
could exceed the classical value by a factor of around 10^. Under 
these circumstances a fluid model can be appropriate even for 
relatively high frequency waves. This may also be true in the up- 
per chromosphere, where the plasma is both cooler and denser 
(and consequently much more collisional) than in the corona. 

It is well known that the electron inertia term in the gener- 
alised Ohm's law becomes comparable to the MHD terms when 
the system lengthscale approaches the electron skin depth, 5e, 
which in a low beta plasma can exceed the ion Larmor radius. 
Moreover for perturbations with frequencies approaching the ion 
cyclotron frequency Q,, the Hall term in the generalised Ohm's 
law becomes as important as the MHD terms when the length- 
scale of the system approaches the ion skin depth 5, (» 6e)- 
When the introduction of the Hall term into Ohm's law is the 
only modification made to the otherwise standard set of MHD 
equations, we may refer to this as a Hall MHD system. 

Hall MHD has been found to be important for a number 
of fundamental plasma processes. For example, in magnetic re- 
connection studies, |Birn et aL] ( |2001| l found that all models 
which included Hall dynamics returned indistinguishable recon- 
nection rates, concluding that the inclusion of the Hall term is 
the minimum requirement for fast reconnection (for a summary 
of Geospace Environmental Modelling (GEM) challenge results, 
see |Bim & Priest| ( [2007] l). 

High frequency waves (i.e. with frequencies oj ~ Q.,), have 
been observed in a range of astrophysical plasma systems, for 



example in the solar corona (summarised in Marsch ( 2006 1) and 
in situ, at the Earth's bow shock ( Sckopke et al.|[T990 l. When 
oscillations in this frequency range are excited in collision- 
dominated plasmas, such that the collisional mean free path is 
less than the wavelength, it is then appropriate to use a Hall 
MHD model. 

The goal of our work is to determine the extent to which 
the main consequences of phase-mixing (wave dissipation and 
plasma heating) are afifected solely by the addition of the Hall 
term to Ohm's law. To do this, we first investigate the damping 
rate of a uniform plasma using Hall MHD (Section |2]i. Phase- 
mixing is then included, by allowing the equilibrium density to 



vary (Section 2.3 1. Numerical simulations of a Hall MHD sys- 
tem, with various density profiles and Hall term strengths, are 
described in Section [3] We interpret these results and present 
conclusions in Section]?] 



2. Wave Damping Analysis 

The Hall MHD form of Ohm's Law is 

E + vxB = 77J-i j X B, 

ne 

where E and B are electric and magnetic fields, v is the single- 
fluid plasma velocity, j is current density, 77 is electrical resis- 
tivity, n is plasma number density and -e is electron charge. 
Linearising the induction equation corresponding to this form 
of Ohm's law, together with the equation of motion (neglecting 
the pressure gradient force), for a plasma with a uniform equilib- 
rium field Bo = Boy for a constant Bq, uniform number density 
and resistivity, and zero equilibrium flow, gives: 

1 



■Q¥J 

— !- = V X (V X Bo) + ^V^Bi 

at Ho fione 



'di 



1 



Pojuo 



(VxBi)xBo, 



Vx[(VxBi)xBo](l) 



(2) 



where Bi [= {Bx, 0, B-)] now represents a transverse perturbation 
to the equilibrium field, po is equilibrium mass density, and ^q is 
the permeability of free space. Note that although, for simplicity, 
we neglect plasma pressure, by assuming y6 = in this analyt- 
ical treatment (as Alfven, whistler and ion cyclotron waves are 
all incompressible in a linear regime), the numerical simulations 
presented in Section[3]incorporate plasma pressure, i.e. /? ^ 0. 
By inserting (J2]i into the time derivative of (fT) in the usual man- 
ner, the linearised Alfven wave equation is then modified: 



d^B, 



= ca 



,d^B, 

5y2 



fio 



- iCA^i 



dB, 

dy^\~dt j' 



(3) 



where we have expressed transverse field perturbations in the 
form of a complex variable Be - B, + iB^, and introduced 
the Alfven speed, c^ - Bq/ yfpofiQ, and the ion skin depth, 
6i - cjiDpi - c ■\jmi£o/ne^, where nij is the ion mass, c is the 
speed of light, and eo is the permittivity of free space. 
Seeking wave-like solutions of the form exp (/ [ky - a>t]) allows 
us to form a dispersion relation to express perturbation frequen- 
cies, (L), as a function of wavenumber k. By considering only 
the regime of weak damping {rfk^ I jj^ca^ <k 1), we obtain two 
separate solutions depending on the size of the parameter k^6^. 
Taking first the case ofk^5^ <K 1, we make use of a simple Taylor 
expansion to find (to leading order in k^5j): 

k^ irjk^ 

(X) = +CAk ± CA5i— - - — , 

2 2//0 

which describes a shear Alfven wave, modified by Hall efifects 
and subject to resistive damping. 
Considering the opposite case, when k^5^^ » 1, we find: 



(4) 



irjk 



^^+c S- + C 6- 1 + ^^ + — 

' 2 ~ '^ ' 2 V //oCa5; k^6^ 2yL(o 



Including the next term in the Taylor expansion, we find two 
distinct forward propagating solutions: 

Ca it] 



0)+ 



^i /"o5. 



(5a) 
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2 ^^^^ 
CO- - CA^ik , 



(5b) 



where in this limit, we now obtain a combination of whistler 
and ion cyclotron (i.e.) waves, both subject to a form of resistive 
damping. 

2.1. Long Wavelength Hall MHD Regime (Uniform) 

We can examine the effect of this difference in behaviour of 
both k^6J regimes, by focussing on the evolution of an initially 
Gaussian pulse (of width cr, and amplitude B,), which is allowed 
to travel along the equilibrium field, taking the form; 



B,(y, 0) = B^ + iB, = B,- exp I - ^ j. 



(6) 



Our complex variable can be interpreted as a Fourier integral, 
evolving in time as: 



f{k)exp{i{ky-ajt)]dk, 

oo 

with f(k) determined by the initial conditions: 



(7) 



1 r°° B (T I k^rr^ ^ 
f(k) = — B,(y, 0) exp (iky) dy = ^— exp — |, (8) 



where we have used the standard result ( Abramowitz & Stegun 



1964 Eq. 7.4.6): 



Xoo 
exp l-lix^\ cos (bx) dx 



1 [^ I b^ 

2 VyS^^^r 4^ 



(9) 



We also use this result to evaluate the integral in Eq. ^ in the 
limit k^6^ <K 1, finding: 



Bc-Yu 



I exp 



-(y + CAtf 



(10) 



where the summation is over forward- and backward 



propagating waves. Eq. ( 10 1 describes a pair of pulses travelling 



in opposite direction at approximately the Alfven speed, which 
are damped by finite resistivity and circularly polarised. We can 
also calculate the contribution to the total energy per unit area in 
(x,z), £>Btot^ made by the magnetic energy per unit area in (x,z) 
of both pulses (£b ) as follows: 

&BroT = :^ (Bo + Bi)2 flfy = -i^ + £b, . 

Since the magnetic perturbation is transverse to the equilibrium 
field. Bo ■ Bi = and hence 

1 r 1 r°° 

£b, = ^ ^i^dy = — B,B*Jy. (11) 

2;U0 J 2llQ J_^ 

Many of the factors in BcB* will cancel upon integration, hence 
the energy associated with the pulse (£s^) evolves as: 






V^o-B^ 



4fio(l + r]t/jj.QO-^) 



.2^l/2 



1 + exp - 



ca t 



0-2 + T^t/jJo 



■ (12) 



Thus after a short initial transient phase (essentially the time 
taken for an Alfven wave to travel a distance equal to the ini- 
tial pulse width, cr), we recover a power law decay (oc r'''^ for 
t » fj.Qcr^/ri) in the energy associated with the pulse. This ex- 
pression (Eq.fTS]) is compared with several numerical simulation 
results in Section |3] and can be seen in Fig. |2] It is straightfor- 
ward to show that the expression given by Eq. ( [T2| for the Hall 
MHD long wavelength (k^dj ^1) regime is identical to that 
found for an initially Gaussian pulse in the MHD limit. 

2.2. Short Wavelength Hall MHD Regime (Uniform) 

Turning to the opposite limit, k^6f » 1, the perturbation fre- 
quencies (J5]l comprise of a combination of resistively damped 
whistler and ion cyclotron waves. We can again describe the 
pulse evolution associated with each separate wave branch in this 
limit, in the manner described previously (Section [2T| , again for 
an initially Gaussian perturbation. Beginning with dTli, and with 
the same initial conditions (J8]l as the previous limit, we find that 
the whistler wave calculation proceeds similarly to that of the 
previous section, however the i.e. wave (being independent of 
wavenumber) differs somewhat: 



B': 



B" 



Bicr 

y/2^ 



e''y 

Jo 



exp 



.2 7,2 



o-^k 



rjk^t ,2 1 „ 

icfydik t\dk, 

juo 



(13a) 



B;Cr ( iCAt Tit 

: exp' 



2lT 



m^f) 



/-*CC 

Jo 



exp I — ;; — \dk. (13b) 



Evaluating the integrals in Eq. ( 13 1, using Eq. (J9]), we obtain: 



B"' = 



'"■ = :exp p^ -f 



(14a) 



B'^' = Bi exp 



ICAt 

'~6~ 



It 



exp 



y 

'2o-2 



(14b) 



In this short wavelength (k^6^ » 1) regime, the peak of the 
pulse now no longer propagates, but decreases in amplitude. The 
right circularly polarised component of the pulse rapidly broad- 
ens, due to the high whistler speed, and damps algebraically at a 
rate similar to that found in both the MHD and long wavelength 
(k^S^ <^ I) Hall MHD regimes. The left circularly polarised (ion 
cyclotron wave) component, on the other hand, damps exponen- 
tially. It should be noted that this damping arises from resistive 
dissipation, and as such should be distinguished from the kinetic 
ion cyclotron damping arising from wave-particle interactions. 
We may, again, calculate the energy associated with each 



solution ( 14 1, using (111, where we still only obtain transverse 
perturbations, and hence Bi ■ Bo still makes no contribution to 
the energy. In this limit, we obtain an expression for the energy 
associated with the individual whistler and i.e. wave branches: 






B^o-y/^ 

2yUo 



1 + 



2r]t 



-1/2 



(15a) 
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ei- 






2)Uo 



(15b) 



Thus, for waves in the k^6] » 1 regime, we no longer see 



the initial transient phase seen previously in (12 1, and for long 



timescales the algebraically-damped whistler contribution to the 
wave energy is dominant over the exponentially-damped contri- 
bution from the ion cyclotron wave. 

2.3. Wave Damping and Phase-Mixing in a Non-Uniform 
Plasma 

By now allowing the equilibrium plasma density to vary in a 
direction perpendicular to both the direction of the equilibrium 
field (y) and the direction of initial perturbation (z), we can in- 
vestigate what effect the Hall term has on the dissipation rate 
in a non-uniform plasma. When the gradients in the jc-direction 
are sufliciently large, and the effects of viscosity are negligible, 
the linearised Alfven wave equation in the MHD limit takes the 
form | |Hoodetal.| ( |2002l l]: 



5^ 



2 d^B, rj d^ IdB- 



(16) 



The variation in Alfven speed ca{x) - Bq/ ^j|lQp{x) causes 
steep gradients to build up in the direction of the inhomogene- 
ity which, in turn, significantly enhances resistive damping in 
the regions where the inhomogeneity is greatest. |Hood et al. 
( 2002| l used a multiple time-scale analysis to derive from Eq. 
( [T6) l a one-dimensional diffiision equation whose solutions can 
be expressed in terms of the Alfven speed gradient ca'{x) - 
dcA{x)/dx. For the case of the initially Gaussian pulse defined 
by Eq. (J6]l, the forward-propagating solution takes the form: 



B, = 



Bi 



2^Jl+CA'^T]ty3flQO-' 



exp - 



(y - CAtf 



2 [cr2 -H Ca' VV3a(o] 



(17) 



We can evaluate the perturbed magnetic field energy per unit 
length in the z-direction Sg for this case by integrating B^/2jUo 
over a finite distance xq < x < xi in the x-direction and from 
minus to plus infinity in the y-direction, taking into account 
the presence of both forward- and backward-propagating pulses. 
The y-integration can be performed analytically, yielding: 

B^o-yfJi r" dx 



oH _ 



8/^0 J,-„ (1 



B]cr^lli 
~8^ 



(l+c^'2;7fV3/zo^2)i/2 

exp(-^2+,y2^,3/3^J 



^0 Ja-o (I+Ca'^TII 



dx. 



(18) 



The integrals in this expression can be readily evaluated numeri- 
cally, to allow comparison with our numerical simulations in the 
next Section. 

3. Simulation Results 

The system was modelled numerically using a two dimensional 
version of a Lagrangian remap scheme (LareXd), described by 



|Arber et al.| ( |2001[ ), which includes an optional Hall physics 
package to incorporate the Hall term into the standard MHD sys- 
tem of equations, seen here in normalised dimensionless form: 



dp 

dt 



4- V ■ (pv) = 0, 



p( -^ -H (v ■ V) v) = (V X B) xB - V/?, 



m 

dt 



V X (v X B) - V X (77V X B) - AjV x 



(V X B X B) 



P\^t^^'^-^'^^\^-P'^-^"''^j'' 



for dimensionless mass density p, pressure p, magnetic field 
strength B, fluid velocity v, internal energy density e, resistivity 
T] (the reciprocal of the Lundquist number) and ion skin depth 
U^ 6ilh). 

In Section [Tl two particular regimes were identified where 
collisional models, such as this, may be more appropriate to 
describe the plasma behaviour than collisionless treatments. 
Normalising temperatures {Tq) and densities (no) using flaring 
coronal values (ro ~ 2 x lO^K and n^) x lO'^m"^) or upper 
chromospheric values (Tq x 2 x lO'^K and no ~ 2 x lO'^m"-') 
places the simulations firmly within these regimes. Specifying a 
low plasma beta (J3 = 0.01) in the simulations fixes the normalis- 
ing magnetic field strengths to Bo « 1 18G in the flaring corona, 
or Bo ~ 17G in the upper chromosphere. This also determines 
the effective size of several fundamental plasma parameters, out- 
lined for these normalising values in Table [T] 

Table 1. Approximate Lengthscales 



Parameter 


Flaring Corona 


Chromosphere 


electron gyro-radius (r<,) 


0.26cm 


0.18cm 


ion gyro-radius (r,) 


11cm 


7.9cm 


electron skin depth (Si,) 


5.3cm 


3.8cm 


ion skin depth ((5,) 


2.3m 


1.6m 


classical mean free path (/1,,,/p) 


30km 


1.5m 



Notes. Approximate values of fundamental plasma parameters, cal- 
culated using normalisation values. In the flaring corona, these were 
density ii ~ 1 x lO'^m"-', temperature T ~ 2 x 10*K and mag- 
netic field strength B ~ 118G, while in the chromosphere, values of 
« ~ 2 X 10'^m-^ r ~ 2 x lO^K and B ~ 17G were used. 



The parameter /i, (which controls the effect of the Hall term 
in our simulations) was initially set equal to 0.0072. Using this 
value of Aj, together with the ion skin depths listed in Table 
[T] implies a normalising lengthscale Iq = 0.3km in the flaring 
corona, and Iq - 0.2km in the chromosphere. Given that the 
mean square wavenumber for a Gaussian pulse of width cr can 
be found using {k^} = 1/cr^, a choice of width cr = 0.1 places 
the simulations firmly within the long wavelength regime (dis- 
cussed in Section : 



2.1 1, as {k'^}6f =s 0.005 <; 1. By relating sim- 



ulated perturbation frequencies (a)) to the perturbation size (by 
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assuming u ~ kc^), our choices for /I, and cr perturb the simu- 
lations with frequencies which have begun to approach the ion 
cyclotron frequency (w/Q,- ~ kdi ~ Aj/cr ~ 0.07). 

As pointed out by [Craig & Litvinenko (20021, the resistiv- 
ity could be as much as a factor of 10^ higher than the classical 
Spitzer value in the flaring corona. Using this enhancement fac- 
tor and the normalisation described above for flaring conditions, 
we find J] = 0.0005. Note that using the same value for rj and 
the chromospheric normalisation, the enhancement of the resis- 
tivity as compared to the Spitzer value, is less than 10^. An en- 
hancement in the resistivity implies a corresponding reduction of 
the mean free paths quoted in Table [T] These reduced mean free 
paths are then much smaller than the typical wavelengths of the 
modes under consideration, so that Hall MHD is an appropriate 
model to use. 

In order to study phase-mixing, we allow the equilibrium 
density to vary along x, with the form: 



Table 2. Density Enhancement Values 



p(x) 



1 



(19) 



(1 -a' + Q'COs(2;rjc))'^^' 

chosen for constant density at the edges, p(0) - p{l) - 1, and a 
central increase in density controlled by a steepness parameter, 
a. We also vary the specific internal energy density of the sys- 
tem, e, allowing us to define a constant plasma pressure in terms 
of the plasma beta, /?(= 0.01), and the magnetic pressure, as: 



e(x) 



PW 



P ^p(xMx)(y- 1) 



m^ 



2p(x)(r-i)' 

We set up a constant equilibrium field to reflect the analytical 
setup, with Bo(= [Bq^, Boy, Bq^]) - [0, 1,0], which we perturb 
in the form outlined by Eq. (|6]l. As our analytical work is based 
on a linearised system of equations, we fix the pulse amplitude, 
Bi{- 0.0005) to be small enough that we minimise non-linear 
eflfects. A 2000 x 2000 grid was used to fully resolve the effects 
of a density gradient in x and potential whistler/i.c. effects in y. 
As the density and energy density profiles are symmetric, peri- 
odic boundary conditions are used in x, whilst the y boundaries 
are left open, and the simulation terminated before the pulse 
reaches the boundaries. The range of the numerical box (set, for 
the k^6j <K 1 simulations, at -10 < j < 10, < x < 1) has a 
large range in y, in anticipation of the pulse propagation. 

Our initial investigation sought to recover, using Lare2d, 
the behaviour of the pulse and its evolution along the equilibrium 
field, in the long wavelength Hall MHD regime (expressed in Eq. 
10 1. For this investigation, selecting a steepness factor of a = 



in our equilibrium density expression (Eq. 19 1 provided the re 



quired uniform equilibrium density and uniform internal energy 
in the simulations. We found that the evolution of the pulse in 
the numerical scheme exactly matched the equivalent analyti- 



cal expression (Eq. 10 1, both in the MHD and long-wavelength 
Hall MHD (Ai - 0.0072) limits. We also sought to recover the 
expression for the evolution of perturbed magnetic energy asso- 
ciated with the pulse in the long wavelength Hall MHD regime 
(Eq. [T2]i numerically. This too showed excellent agreement with 



Steepness Parameter (a) 


Pinax Pmiii 


1/10 


0.12 


9/50 


0.25 


5/18 


0.50 


3/8 


1.00 


4/9 


2.00 


15/32 


3.00 



Notes. Values of a are chosen to provide specific sizes of density en- 
hancement, p„oj- - p„„„, at the centre (x = 0.5) using Eq. jl9[). 



the derived expression, and a comparison of the numerical and 
analytical energies can be seen in Fig. [2] 

Moving to the non-uniform equilibrium density simulations, 
we then selected a range of specific equilibrium density enhance- 



ments through the steepness parameter a according to Eq. ( 19 1 



The specific values of a chosen, and the corresponding density 
enhancements are listed in Table |2] 
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Fig. 1. Comparison of predicted amplitude fall-off' for an ana- 



lytical treatment of MHD phase-mixing (Eq 17 - solid curve) 



with tracked pulse peak of Lare2d simulations. Displayed are 
the results of the MHD (crosses) and long wavelength Hall MHD 
(diamonds) simulations for three different density gradients (for 
corresponding density enhancements, see Table|2]i. 
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As discussed previously, |Hood et al.| ( |20(J2) l obtained an ex- 
pression describing the evolution of a pulse as it undergoes phase 



mixing in the MHD limit [Eq.( 17 1]. We sought to test the valid- 
ity of this expression for a range of equilibrium density gradi- 
ents using Lare2d, by tracking the height and location of the 
peak of the pulse as it travelled along the equilibrium field, and 
comparing this with the pulse amplitude damping rate found by 
|Hood et alT] ( |2002| l. We found that the MHD simulations with 
large equilibrium density increases displayed good agreement 
with the proposed damping rate, but this agreement was lessened 
by reducing the equilibrium density enhancement. Furthermore, 
we repeated the simulations for the long wavelength Hall MHD 
regime, again finding that the amplitude evolution of the pulse 



was in agreement with the expression proposed by Hood et al 



( |2002| l, but only in the cases where the density enhancement was 
sufficiently large. These results are summarised in Fig.[T] 

We also evaluated numerically the integrals in Eq. (flSll to 
investigate how the magnetic energy associated with a pulse, 
which undergoes MHD phase-mixing, evolves in time (see 
Section [23] l. As we have shown, the expression of |Hood et al.| 
( 2002| l displays good agreement with the simulations in the cases 
of steepest density gradient. We therefore selected the steepest 
of our density gradient cases (a = 15/32), prescribing a par- 



ticular equilibrium density p(x) (from Eq. 19 1, and hence ob- 
tained an expression for the gradient in Alfven speed, c/{x), re- 



quired for the numerical integration of Eq. ( 18 i. Using Lare2d, 
we then compared this with the evolution of magnetic energy of 
an identical pulse and equilibrium density gradient, also using 
MHD. This comparison can be seen in Fig. [2] together with the 
magnetic energy evolution of an identical pulse for the uniform 
density case, and the analytical expression derived for the uni- 
form case, both in MHD. In the uniform density case there is 
a very close agreement between the analytical and Lare2d re- 
sults. When a density gradient is present, the field perturbation 
energy starts to decay slightly earlier in the Lare2d simulations 
than it does in the analytical result. In all cases, there is a short 
transient phase after f = in which the field perturbation energy 
drops to half its initial value and the flow energy reaches approx- 
imate equipartition with the field energy. As discussed in Section 
|2j the duration of this transit phase is essentially determined by 
the Alfven propagation time across the initial pulse width. 

Having studied the behaviour of the perturbed magnetic en- 
ergy evolution in the uniform and steep density gradient lim- 
its using MHD simulations, we then investigated the energy re- 
sponse for a variety of different density gradients, and at different 
skin depths. For clarity, we also included the perturbed internal 
energy evolution of the pulse. Our first set of results compare 
the evolution of the pulse in MHD, with that of the Ai - 0.0072 
simulations, corresponding to the long wavelength Hall MHD 



regime, discussed in Section 2.1 These results are seen in Fig. 

m 

FigurelSlclearly displays differences between the MHD and 
Hall MHD cases, for density gradients which are neither uniform 



Lare2d, MHD, uniform 

anolytic, MHD, uniform 

LareZd, MHD, K = 15/32 

Hood(2002) integrated, c< = 15/32 



-m — X — ^-5^_ 



Fig. 2. Temporal evolution of magnetic energy perturbation as- 
sociated with initially Gaussian pulse, computed analytically for 
cases with uniform equilibrium density using Eq. (fTTli (stars). 



and strong equilibrium density gradient using Eq. ( 1 8 1 (dotted 
curve). The solid and dashed curves show the corresponding nu- 
merical results, obtained using Lare2d. 



nor very steep. To gain insight into the reasons for these differ- 
ences, we have also tracked the evolution of the pulse profile in 
the jc-direction as it travels at the Alfven speed of the steepest 
density region. In the Hall MHD cases, we observe consistently 
smaller amplitude gradients. As an example of this, Fig.|4]shows 
snapshots of MHD and Hall MHD pulse profiles in simulations 
with a = 5/18. One might have expected the dispersive effects 
of the Hall term to induce a cascade into short wavelengths, 
which could in principle invalidate our neglect of kinetic effects. 
However, Fig. [4] shows that the opposite effect occurs: the Hall 
term actually reduces the gradients resulting from phase mixing. 

We have demonstrated analytically that the plasma response 
dramatically differs in the short wavelength (k^6f » 1) Hall 
MHD limit, from that of the MHD and long wavelength Hall 
MHD limits (Section l2|. Our final investigation sought to be- 
gin to bring out the behaviour of this limit by increasing A/ by a 
factor of 10. In doing so, we also reduce the associated enhance- 
ment factor in resistivity required to justify setting rj - 0.0005 in 
the code. The enhancement factor is reduced to 10^ for the flar- 
ing corona, whilst for the upper chromosphere the effective re- 
sistivity is less than an order of magnitude larger than the equiv- 
alent Spitzer value, using the values outlined earlier in Section 
|3] In these simulations, we also increased the size of the nu- 
merical box (to -30 < y < 30, < x < 1, accounting for 
larger dispersive effects and the faster whistler wave component) 
and the number of gridpoints (now at 6000 x 2000, maintaining 
the same resolution). The perturbed frequencies in these simu- 
lations are now a large fraction of the ion cyclotron frequency 
(dj/Q., ~ Aj/cr ~ 0.7). Fig. [5] compares the response of the per- 
turbed magnetic and internal energies in simulations of a Hall 
MHD plasma with /I, = 0.072, with that of an MHD plasma, for 
three density gradient cases which are neither uniform, nor dra- 
matically varying. It is clear that in this Hall MHD regime, there 
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Fig. 3. Evolution of field perturbation energy in MHD (solid curve) and A, = 0.0072 Hall MHD (dashed curve) simulations for six 
values of density steepness parameter a. The internal energy is also plotted for the MHD (dot-dash curve) and /I, - 0.0072 Hall 
MHD (dotted curve) simulations. 
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Fig. 4. Comparison of two snapshots of slices (in x) through the 
pulse amplitude, at location (in y) of maximum phase-mixing. 
Snapshots shown are for long wavelength Hall MHD (solid) and 
MHD (dashed) cases, taken at the same time, for an identical 
initial pulse and with density steepness a = 5/18. 



is a strong reduction in the damping rate compared to the MHD 
limit. 



4. Discussion and Conclusions 

The behaviour demonstrated by the analysis of a uniform plasma 



in the long wavelength Hall MHD regime (Section 2. 1 1 was fully 
recovered by the Lare2d simulations. Both the simulated pulse 
evolution and its associated magnetic energy are indistinguish- 



able from the expressions derived analytically (Eq. ( 10 1 and Eq. 
( [TT| i respectively). The simulations show that the MHD results 
remain linearly polarised, whilst the long wavelength Hall MHD 
pulse becomes circularly polarised. However, both the simula- 
tions and the analytical results tell us that this has no impact on 
the magnetic energy associated with the pulse, which damps at 
the same rate in both regimes. 

Considering now the results of the investigation concerning 
a non-uniform equilibrium density. Fig. [T] shows that the am- 



plitude damping rate described by Hood et al. (2002i is well 
matched to the rate of damping of the peak of an MHD pulse 
when the gradient in Alfven speed is sufficiently steep. Whilst 
not a particularly surprising result, as the strong phase mixing 



limit required for Eq. ( 16 1 to be valid applies in this case, it is 



somewhat more surprising to find for sufficiently steep Alfven 
speed gradients, that the long wavelength Hall MHD regime also 
conforms to this amplitude damping rate. For shallower gradi- 
ents in Alfven speed, Fig.fllalso demonstrates that this damping 
rate quickly becomes inappropriate for the k^6j <K 1 Hall MHD 
simulations, whilst giving better agreement with the MHD evo- 
lution of the peak of the pulse. This is partly due to the fact that in 
the Hall regime, the pulse is no longer linearly polarised, quickly 
losing amplitude to the other field components. The Hall regime 
is also subject to significant dispersive effects, which reduce the 
amplitude of the pulse as it spreads along the equilibrium field. 
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Fig. 5. Evolution of field perturbation energy in MHD (solid curve) and /I,- - 0.072 Hall MHD (stars) simulations for three interme- 
diate values of density gradient parameter a. Also plotted is the internal energy for the MHD (dot-dash curve) and At - 0.072 Hall 
MHD (crosses) simulations. 



However, for sufficiently steep density gradients, phase-mixing 
becomes the dominant cause of amplitude dissipation. 

By comparing the energy evolution associated with our nu- 
merical simulations of an MHD pulse to that found by numeri- 
cally integrating an expression for the pulse evolution given by 
|Hood et aL ( 2002| l in Fig.l2] we see that, apart from a slight initial 
delay, the two results appear closely matched. The reason for the 
difference between the two is simply due to the treatments used. 



The linearised Alfven wave equation ( 16 1 used in the treatments 



of |Heyvaerts & Priest| ( |1983| l and |Hood et al.| ( |2()02| i retains only 
the x-derivatives, on the basis that derivatives along the field are 
small compared to those across it. However, in our simulations, 
the (initially uniform) pulse propagates a finite distance in which 
the y-derivative causes resistive damping along the equilibrium 
field. This continues until the gradients in x build up sufficiently 
to allow phase-mixing to become the dominant damping mech- 
anism. The brief initial period of longitudinal damping is not 



included in the treatment of Hood et al. ( 2002) l, resulting in the 
discrepancies in behaviour seen between the energy evolution of 
the two non-uniform scenarios shown in Fig. [2] 

Our final investigation sought to compare the response of 
the magnetic and internal energies associated with the pulse 
for different equilibrium density gradients and at different skin- 
depths. Fig. [3] shows that for shallow density gradients, both 
MHD and long wavelength Hall MHD simulations display sim- 
ilar behaviour to that seen in the uniform case. For the steep- 
est density gradients considered, again the long wavelength Hall 
MHD and MHD simulations exhibit similar behaviour, suggest- 
ing that MHD phase-mixing dominates the energy evolution (see 
Fig. [2]). However, we see significant departures from the MHD 
behaviour of both perturbed magnetic and internal energies asso- 
ciated with the long wavelength Hall MHD regime, in the cases 
where the density is neither uniform nor sharply varying. In these 
cases, the inclusion of the Hall term causes the damping rate of 
perturbed magnetic energy to be reduced. Consequently there is 
a slower increase in the internal energy of the pulse than that seen 



in MHD: phase-mixing of the plasma causes less rapid plasma 
heating in the Hall MHD regime than in the MHDregime. 

Evidence to support this can be seen in Fig. HI where we 
plot a slice (in x) through pulse amplitude at the location (along 
y) of maximum phase-mixing, both for the Aj - 0.0072 and 
MHD simulations. The figure clearly demonstrates that in the 
Hall case, the gradients in x are smaller than those in the MHD 
case. The whistler component of the Hall MHD pulse displays 
highly dispersive behaviour, which increases with skin depth. 
This dispersion spreads the pulse envelope along the equilibrium 
field as it travels, introducing higher wavenumbers into the pulse 
but only in the equilibrium field direction. In the density gradi- 
ent direction, the group dispersion of the wave has the effect of 
reducing the damping rate by making amplitude gradients across 
the field (in x) smaller, reducing the efficiency of phase-mixing 
and its damping of the magnetic energy of the pulse. 

This behaviour is also seen in the higher skin depth Hall 
MHD simulations. In the Aj - 0.072 simulation results, seen in 
Fig. |5] there is now a very significant reduction in the damping 
rate of the magnetic energy, even in the steep density gradient 
case, which, in the long wavelength Hall MHD regime has begun 
to converge to the MHD results. As the skin depth is increased, 
the equilibrium density gradient must also be increased in order 
to recover a rate of energy dissipation comparable to those seen 
in the long wavelength Hall MHD and MHD hmits. 



In contrast to the results presented here, Tsiklauri et al. 
C2005') and 'Bi an & Kontar| ( |2010 ) found that by extending the 
MHD treatment of phase mixing to include kinetic effects it was 
possible to demonstrate an enhanced rate of wave dissipation 
rather than a reduced rate. Specifically, Tsiklauri et al. ( 2005| l 
observed Alfven wave damping in a collisionless particle-in- 
cell simulation of a plasma with an equilibrium density gra- 
dient, the amplitude decay law being similar to that found by 
'Heyvaerts & Pries^t] ( |1983| l using an MHD treatment. ' Bian &| 
,Kontar, ( ,2010) used drift-kinetic theory (applicable to mode fre- 
quencies CO « Q.i) to argue that this result can be attributed to 
the generation of a parallel electric field E\\ associated with the 
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presence of finite perpendicular wavenumbers and consequent 
mode conversion to kinetic Alfven waves; because /sy is finite 
the wave energy can be dissipated via Landau damping. The re- 
sults of our fluid study show that extensions of phase mixing 
theory to include non-MHD effects do not necessarily lead to an 
enhanced damping rate in the case of waves with u) ~ O, and in 
the presence of collisions. In view of the differences between the 
physics assumptions of Hall MHD on the one hand and those of 
collisionless Vlasov or drift-kinetic theory on the other, it is not 
particularly surprising that the two models produce different re- 
sults with regard to phase mixing. The collisionless regime is the 
relevant one for high frequency waves in the corona when the re- 
sistivity is close to classical. Our calculation, on the other hand, 
may be more applicable to either flaring conditions in the corona 
or the upper chromosphere. In our current study, we have only 
included dissipation through resistivity, neglecting other dissipa- 
tive mechanisms. ^Craig et al. ( 2005| l noted that classical viscos- 
ity could play a significant role in flare energy release. Hence, a 
possible extension of the work reported here would be to exam- 
ine the effects of the Hall term on phase mixing when viscous 
terms are added to the momentum and energy equations. 

In summary, we have described analytically the propaga- 
tion of an initially Gaussian field perturbation along a uniform 
equilibrium field in the presence of resistivity, and the evolution 
of magnetic energy associated with this perturbation. While the 
evolution of the perturbation differs in the MHD, long wave- 
length and short wavelength Hall MHD regimes, the energy 
evolution is the same in the MHD and long wavelength Hall 
MHD cases. In a non-uniform equilibrium plasma, our simula- 
tions show that the damping rate for the energy associated with 
the pulse in Hall MHD is significantly reduced when the Alfven 
speed variation is neither uniform nor sharply varying, compared 
to that of an MHD treatment as the Hall term actually reduces 
the gradients resulting from phase-mixing. Moreover, as the ion 
skin depth is increased, the density gradient needed for MHD 
phase-mixing to dominate the evolution of the pulse (the "strong 
phase-mixing limit") must also increase. 
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